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We discuss the dynamic properties of the square-lattice spin-| XY model obtained 
using the two-dimensional Jordan- Wigner fermionization approach. We argue the rele- 
vancy of the fermionic picture for interpreting the neutron scattering measurements in 
the two-dimensional frustrated quantum magnet CS2CUCI4. 
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Dynamic properties of two-dimensional quantum spin models have attracted 
considerable attention from both experimentalists and theoretists. Among recent 
experimental studies in this field the neutron scattering on CS2CUCI4 03 E] 
seems to be very intriguing. CS2CUCI4 is a two-dimensional low-exchange quan- 
tum magnet. The two-dimensional character arises owing to the layered crystal 
structure, whereas the in-plane superexchange route mediated by two nonmagnetic 
Cl~ ions provides a small exchange interaction (~ 1 — 4 K) between Cu 2+ ions 
carrying a spin of ^. In each layer the exchange paths form a triangular lattice 
with nonequivalent interactions along chains (J) and along zig-zag bonds (J'). The 
minimal Hamiltonian which determines the magnetic order is the two-dimensional 
Heisenberg Hamiltonian 

H = ^ JSi ■ Si' + ^ J' s i ■ s j; (!) 

where the first term describes the noninteracting chains (J = 0.374(5) meV) and 
the second one describes the interchain coupling (J' — 0.34(3)J). The interlayer 
coupling is small J" = 0.045(5)J and it stabilizes the long-range magnetic order 
below Tn = 0.62(1) K. The neutron scattering measurements in the spin liquid 
phase (i.e. above Tn but below J, J' when the two-dimensional magnetic layers are 
decoupled) clearly indicate that the dynamic correlations are dominated by highly 
dispersive excitation continua which is a characteristic signature of fractionalization 
of spin-1 spin waves into pairs of deconfined spin-i spinons. Linear spin-wave theory 
including one- and two-magnon processes cannot describe the continuum scattering. 
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Two theories of the dynamical correlations in the spin liquid phase in CS2CUCI4 
have so far been proposed. A quasi-one-dimensional approach J5_ starts from an- 
tiferromagnetic Heisenberg chains along the strongest-exchange direction in the 
two-dimensional triangular planes (thus immediately introducing spinon language) 
and perturbatively treats all other terms in the Hamiltonian. Note, however, that 
the coupling between chains is relatively large. Another approach is based on the 
explicitly two-dimensional resonating- valence-bond picture in which spinons appear 
to be two spin-i ends of a broken singlet bond that separate away through bond 
rearrangement jBJ |7| . 

In what follows we discuss the two-dimensional Jordan- Wigner fermionization 
approach to dynamic correlations [Hj which eventually may capture the physics in 
spin liquid phase as well. To be specific, we consider a simpler square-lattice spin-^ 
XY model defined by the Hamiltonian 

OO OG 

H = ^ ^2 [J (sijsf +lj j + sf ,jS^ +1 j) + Jj_ (sfj-sfj+i + s i,j s i,j+i)] ■ (2) 

Here J (Jj_) is the exchange interaction in the horizontal (vertical) direction. We 
are interested in the zz dynamic structure factor 

S zz (k,us) = '^ / '^2exp[i(k x p + k y r)] / dtexp(iwi) 

p=0r=0 J -co 

» • (3) 

Performing the transformation to spinless fermions 9 , 

oo oo 

s i,j = ex P ,j ) d i J > «i J = E 3 ln [/ - * + 1 (.9 - j)] d ), 9 d f,g i ( 4 ) 

/=0(^i) g=0&j) 

and adopting the mean-field treatment of the phase factors which appear in @ we 
arrive at the Hamiltonian 

1 oc oo 

H = 2 EE [^y +3 3 (</'••-< - ^-4fu) + J ^ (4,Ai+i d *Ai + i)] > ( 5 ) 

i=0 3=0 

which can be diagonalized by the Fourier and Bogolyubov transformations. The 
zz correlation function which enters is the density-density correlation function 
in fermionic picture JSJ; it can be calculated with the help of the Wick-Bloch-de 
Dominicis theorem. As a result we get 

S„ (k, u) = J g {C 2 [n q (1 - n q+k ) S (w + A q - A q+k ) 

+ (1 - n q ) n q+k 5 (w - A q + A q+k )] + (l - C 2 ) [n q n q+k <S (w + A q + A q+k ) 

+ (1 - n q ) (1 - n q+k ) 6 (w - Aq - A q+k )]} , (6) 
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Fig. 1. S zz (k, u) of the square-lattice spin-i XY model © with J = 1, J ± = 0.3 at 
zero temperature. Frequency shapes are shown along k x = ~, ky = (a), k x = f , ky = 
(b), fc* = fc H = (c), k x = 7T, fc M = (d), fc x = 0, k y = f (e), = f, fc s = f (f), 

k x = 7T, fc y = | (g), fcs, = 0, fc„ = 7T (h), k x = f , fc y = 7T (i), fca; = 7T, fcy = 7T (j). 



where C = cos 7q +^ 7q , cos 7k = J ^^° sfc " , A k = J J 2 sin 2 + J 2 cos 2 fc y and 

n k is the Fermi function. The zz dynamic structure factor (j^J similarly to the 
one-dimensional cases ^U] is conditioned by particle-hole pair excitations. This 
two-fermion excitation continuum implies nonzero values of S zz (k^ oj) only in the 
restricted k-w region with a sharp high-frequency cutoff. Within the excitation 
continuum S zz (k,u>) exhibits several washed-out excitation branches, in particular, 
the spin wave [BJ. The representative scans throughout the Brillouin zone of the 
zero-temperature S zz (k,uj) plotted in Fig. show an extended structure and, at 
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least at a qualitative level, strongly resemble some features of the frequency shapes 
measured by neutron scattering in CS2CUCI4 (see Fig. 6 of Ref. (J)- 

To summarize, we have computed the zz dynamic structure factor of the square- 
lattice spin-i XY model using the two-dimensional Jordan-Wigner fermionization 
approach. The dynamic structure factor exhibits both sharp particle-like peaks 
and extended continuum profiles in the constant wave vector frequency shapes 
which arise owing to particle-hole excitations. Such features which have been ob- 
served experimentally in CS2CUCI4 appear immediately within the two-dimensional 
Jordan-Wigner fermionization approach that is obviously an attractive feature of 
this scheme. One may also expect that in the limit of small interchain interaction 
this approach works well since it contains the exact result if J±_ — (or J = 0). 
Note, however, that this is not a perturbation theory with respect to the interchain 
interaction. The main shortcoming of the exploited approach is a poorly controlled 
mean-field treatment of the phase factors that arise after inserting (0J into J2J- New 
theoretical work is needed to go beyond this approximation JT] . Another challenge 
is to clarify whether such a theory is capable of capturing the physics in spin liquid 
phase in CS2CUCI4. For this purpose the calculations should be repeated for the 
Heisenberg model on a triangular lattice (probably with the Dzyaloshinskii- 
Moriya interaction [3]). We have left this problem for future studies. 
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